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Abstract 

Quantum computation using qubits made of two component Bose-Einstein condensates (BECs) 
is analysed. The use of BECs allows for an increase of energy scales via bosonic enhancement, 
resulting in gate operations that can be performed at a macroscopically large energy scale. The 
large energy scale of the gate operations results in quantum algorithms that may be executed at 
a time reduced by a factor of N, where N is the number of bosons per qubit. The encoding of 
the qubits allows for no intrinsic penalty on decoherence times. We illustrate the scheme by an 
application to Deutsch's and Graver's algorithms. 
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Recent advances in semiconductor technology has allowed the realization of single qubit 
decoherence rates as small as ~ 3 x 10 -6 js 7 ] in spin quantum dots. Such decoherence 
rates are a result of the ability to optically control the qubit states using bright coherent 
pulses, resulting in ultrafast picosecond gates. This allows for a vastly increased number of 
gates allowable within a limited decoherence time T 2 . One of the key ingredients for the 
success of these ultrafast gates originates from the effect of bosonic enhancement. The Rabi 
frequency between two levels is given by Q = fj,E/h, where /i is the optical transition dipole 
matrix element and E is the electric field amplitude. Since E is typically a macroscopically 
large quantity, in principle it is possible increase the energy scale (and thereby decrease 
the time scale) of the single qubit gate by simply increasing the laser power. This bosonic 
enhancement is ubiquitous in many areas of quantum optics and condensed matter physics, 
in phenomena such as laser physics, resonance flourescence and Bose-Einstein condensation 
(BEC), where quantum coherence is "amplified" when a large number of bosons are present 
in the system. Our aim in this work is to bring quantum computation to a macroscopically 
large energy scale, much in the same way that lasers bring the quantum coherence of light 
to a macroscopic level. 

Although the energy scale for single qubit rotations may be amplified in the way described 
above, for two qubit gates, a straightforward extension of this idea is typically problematic. 
For example, consider a two qubit gate connected via a photonic quantum bus (9), [lo| . 
Introducing a large population of photons into the quantum bus typically results in additional 
sources of decoherence, due to the leakage of the photons out of the system. Thus typically 
very weak excitations are assumed in the quantum bus to minimize such effects. Here we 
describe an alternate way of increasing the energy scale of the interactions between qubits 
by populating each qubit with a large number of identical bosons. The large number of 
bosons on each qubit possesses a bosonic enhancement that increases the energy scale of the 
control Hamiltonian, which equates to faster control operations. 

n 

Consider a bosonic system with two independent degrees of freedom [1], such as two 
hyperfine levels of an atomic Bose-Einstein condensate (BEC) or spin polarization states of 
exciton-polaritons Such two component BECs have been realized and controlled in atom 
chip systems, with the demonstration of spin squeezing {^J. Denote the bosonic annihilation 
operators of the two states as a and b, obeying commutation relations [a, a'] = [b, b*] = 1 
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FIG. 1: Two bosonic qubits mediated by a quantum bus. The quantum bus couples transitions 
between levels b n and c n with energy g. Individual pulses coupling levels a n and c n with frequency 
Q n create an adiabatic passage between levels a n and b n . 



is 



8] . We then encode each qubit in a state 



a ,P)) = -=(at£+ptf) N \Q) 



'Nl 



(1) 



where a and are arbitrary complex numbers satisfying \a\ 2 + |/3| 2 = 1 (double brackets 
are used to denote the bosonic qubit states). We consider the boson number N = a) a + tfb 
to be a conserved number. 

The state \a, (3)) can be manipulated using Schwinger boson (Stokes operators) operators 
S x = cvb + tfa, S y = —ia*b + ib^a, S z = a) a — b^b, which satisfy the usual spin commutation 
relations [Si, Sj] = 2ieijkSk, where is the Levi-Civita antisymmetric tensor. In the spin 
language, (QQ) forms a S = N/2 representation of the SU(2) group. Single qubit rotations can 
be performed in a completely analogous fashion to regular qubits. For example, rotations 
around the z-axis of the Bloch sphere can be peformed by an evolution 
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Similar rotations may be performed around any axis by an application of 



H quh = n ■ S = n x S x + n v S y + n z S z 



(3) 
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where n = (n x ,n y , n z ) is a unit vector. Expectation values of the spin are exactly analogous 
to that of a single spin, taking values 

(S x ) = N{a*(3 + a(3*) 
(S*) = N{-ia*f3 + iaf3*) 

(S z ) = N(\a\ 2 -\P\ 2 ). (4) 
The variance of the spins however diminish relative to the maximum amplitudes 

(((S z ) 2 } - (S z } 2 )/N 2 = A\a(3\ 2 /N (5) 

in accordance to widespread notion that for iV — > oo the spins approach classical variables. 

Two qubit interactions can be performed using a generalization of the methods in Ref. 
[lo| . Consider an interaction Hamiltonian 

H hus = y E F « + + 9 E [ F nP ] + F nP] , (6) 

n=l,2 n=l,2 

where F z = c^c — b^b, F + = c^b, uq is the transition energy, and p is the photon annihilation 
operator (we set H = 1). Assuming a large detuning A = tu — tu ^> gy/N, we may 
adiabatically eliminate the photons from the bus by assuming p^p = and we obtain an 
effective Hamiltonian H hus w ^ (F^F^ + F{ F^\ Now consider a further detuned single 
qubit transition according to H pul = Q(t) X^ n =i 2 [ c li a n + H.cl . After adiabatic elimination 
of level c by assuming c^c n = 0, we obtain 

- ^ (^^ + SIS^) . (7) 

The energy scale of the interaction term is then ~ 0(g 2 QN 2 / A 2 ), since the spin operators 
themselves are of order S z ~ O(N). Therefore, by replacing the qubits by the collective 
states ([1]) we have managed to create an interaction term with a boosted energy scale of 
iV 2 . We note that alternative ways of producing two-qubit interactions are also possible, 
such as by using cold controlled collisions [111 ]. 

The consequence of the boosted energy scale of the interaction can be observed by ex- 
amining explicitly the state evolution of two qubits. Before showing this, let us note here 
that the combination of H quh and H int may be combined to form an arbitrary Hamiltonian 
involving spin operators according to universality arguments [12j. By successive commuta- 
tions an arbitrary product of spin Hamiltonians H = Yln=i(^n) m ^ ma y be produced, where 



M is the total number of qubits, a = x,y, z, and m{n) = 0, 1. Using this fact, we consider 
henceforth the interaction Hamiltonian H mt = S'fS'f instead of (J7J), for the simplicity of the 
analysis. As a simple illustration, let us perform the analogue of the maximally entangling 
operation 

e -***i(| t ) + | |) + | I)) = | + y)\ t) + | - y)\ I), (8) 

where | ±y) = e Tl ^ \ f) + e ±l 5| I). Starting from two unentangled qubits, we may apply H iJlt 
to obtain 

e-*wiJ_ JL»iJ_ JL» = JL y ( N )( N 
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where in the second equality we have introduced eigenstates of the S z operator \k) = 
^-4==L=- 1 0) . For gate times equal to t = ir/AN we obtain the analogous state to (JHJ). 
For example, projecting qubit 2 onto the maximum z eigenstates | A: 2 = 0, N) gives the 
states | e± J^ /4 , eT J^ 4 ) ) ? which is the analogue of a Bell state for the bosonic qubits. 

A visualization of the state (JHJ) is shown in Figure [2b. For each z-eigenstate on qubit 
2, there is a state | e ' (JV ^ )7r/4iV ; e t N ^ ^ qubit 1 represented on the Bloch sphere 

entangled with it. The states on qubit 1 are not orthonormal to each other except for the 
extremal states | , e ^ /4 )), but the z-eigenstate on qubit 2 are of course orthonormal. 
Thus the type of entangled state is a continuous verson of the original qubit sequence 
(jHJ), and has similarities to continuous variable formulations of quantum computing [ijj], 
although the class of states that are used here are quite different. We note here that the 
analogue of the CNOT gate can be produced by further evolving (j9j) with the Hamiltonian 
H = NS( - NSI + N 2 for a time t = tt/AN which gives 

—ink2/N i 

<«>) 

Projecting onto the \k 2 = 0,N) states, qubit 1 is left in the states | ± 4g, -^}) which is a 
CNOT operation with the target qubit in the a;-basis for qubit 1. 
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The crucial point to notice in (Q is that a gate time of t = tt/4N was required to produce 
this entangled state, in comparison to the standard qubit case of t — 7r/4. The reduced gate 
time arises due to the boosted energy scale of the interaction Hamiltonian, which in turn is 
due to bosonic amplification of the interaction. This has similarities with other proposals 
for decreasing cooling times by bosonic final state stimulation |4j . 

Despite the widespread belief that for iV — > oo the spins approach classical variables 
according to ([5]), the entangling operation (Q generates genuine entanglement between the 
bosonic qubits. As a measure of the entanglement, we plot the von Neumann entropy 
E = — Tr(plog 2 p) [yl in Figure [2^,. For the standard qubit case (N = 1), the entropy 
reaches its maximal value at t = 7r/4 in accordance with ([H]). For the bosonic qubit case 
there is an initial sharp rise, corresponding to the improvement in speed of the entangling 
operation, but later saturates to a non-maximal value due to the presence of the binomial 
factors in ([9]) biasing the states towards zero spin values. In Figure [2b we show the amount 
of entanglement present at times t = ir/4N, corresponding to the analogous state to (El)- 
We see that at such times there is approximately the same amount of entanglement as for 
the N = 1 for large N, reinforcing the intuition that the e tS i s 2 7T / 4N g a te gives the 

bosonic analogy to the operation ([8]). 

The purpose of boosting the energy scale of the Hamiltonian is defeated if decoherence 
times also decrease when the number of bosons are increased. Naively one might guess that 
a large number of bosons encoding a single qubit would be more susceptible to decoherence. 
This is however is not necessarily always true. As a simple model of dephasing and particle 
loss, consider the evolution of Lindblad form through the master equation 

f = -r^K,K, P ]], (ii) 

n 

where T is the dephasing rate, and A n is the coupling operator to the environment, taken 
to be S* for dephasing and a n , b n for particle loss. For a multiple qubit system, the infor- 
mation in a general quantum state can be reconstructed by 4 M — 1 expectation values of 
(Ji, Sf , Sf, Sf) <g> (J 2 , SI, Si SI) ■ ■ ■ ® (J M , S X M , S V M , S Z M ) [l6]. For the bosonic system, there 
are in general higher order correlations involving powers of operators beyond order one (e.g. 
(S^) 2 ). However, since we perform the bosonic mapping in order to perform a simulation of 
the qubit system, the presence of such expectation values are of no concern in our case. 
Examining the dephasing of the general correlation (Yl n S% ) where a(n) = I,x,y,z, 




FIG. 2: a The entanglement normalized to the maximum entanglment (-E max = log 2 (iV+l)) between 
two bosonic qubits for the particle numbers as shown, b Entanglement at a time t = tt/4N for 
various boson numbers N. c A schematic representation of the entangled state Q. 

we obtain the evolution equation d(Y\ n S% }/dt = — 4TK(Y\ n S% }, which has a solution 

(n^ (n) >ocexp[-4r^]. (12) 

n 

Here K is a constant which is at most equal to M, depending the particular a(n) chosen. 
This equation does not have any N dependence, and in fact behaves identically to the 
qubit case (N = 1). Similar results hold for particle loss. Physically this difference is 
due to the statistical independence of the dephasing and particle loss processes among the 
bosons. Given a qubit algorithm intended for two-level qubits, how does this translate in 
the bosonic system? We have found that for most applications, the procedure amounts 
to: (i) finding the sequence of Hamiltonians required for the algorithm, (ii) making the 
replacement a% — > NS%, cr"a^ —> S^S^, (iii) Evolving the same sequence of Hamiltonians 
for a reduced time t — > t/N. This approach is reasonable from the point of view that we are 
performing the same algorithm except that a higher representation of SU(2) is being used. 
Let us illustrate this procedure with two well-known quantum algorithms with speedups 
over classical algorithms. 
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FIG. 3: a Schematic energy level structure of the Grover Hamiltonian. Rabi oscillations take place 
between the initial x-eigenstate \X) and the solution state \ANS). b Rabi oscillations executed 
by the Grover Hamiltonian for M = 2 for various boson numbers as shown. 



Deutsch's algorithm |l3j. We reformulate the standard qubit version (N = 1) of the 
algorithm in the following form convenient for our purposes. The oracle performing the 
function \x)\y) — > \x)\f(x) © y) is assumed to be one of the four Hamiltonians Hp = 
{0, 2crf , o\a\ + erf — 1, —erf erf + erf — 1} and evolved for a time t — 7r/4, which correspond to 
the functions f(x) = {(0, 0), (1, 1), (0, 1), (1, 0)} respectively. The initial state is assumed to 
be the state (t + 4) t> and a measurement of qubit 1 in the x-basis distinguishes between 
constant and balanced functions via the results (f + I) and (t — i) respectively. 

This can be translated into the corresponding algorithm for bosonic qubits according to 
the following procedure. The oracle is assumed to be one of the following Hamiltonians 
H D = {0, 2NSI , SfS% + NS$ - N 2 , -SfSS + NS% - N 2 }, and we prepare the initial state as 
0)). After evolving the Hamiltonians for a time t = tt/AN, we obtain (up to an 
overall phase) 

e-^ W4 ^|^,^ f ))|l,0)) = |-^,±-J=))|l,0)), (13) 

where + is obtained for the constant cases and — for the balanced cases. A measurement 
of qubit 1 distinguishes the constant and balanced cases with one evaluation of the oracle. 
Grover' s alqorithm. We use the continuous time formulation of the Grover search al- 

n 

gorithm (see sec. 6.2 of Ref. |13|). For the standard qubit case (N = 1), a Hamiltonian 
E G = \X)(X\ + \ANS)(ANS\ is applied to an initial state \X). Here, \X) is the o x n = 1 eigen- 
state of all the qubits and \ANS) is the solution state. Under this Hamiltonian evolution, 
the system executes Rabi oscillations between \X) and \ANS) with a period of t = ir\ / 2 AI 
where M is the number of qubits. 
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The bosonic version of the algorithm can be constructed by first writing the projection 
operators in H as spin operators and following the prescription as described above: 



H G = N 2 f\- 1 + S +N 2 f\- 
J-J- 2 N 11 2 

n=l 



2 

n=l 



1 + f 



(14) 



where we assumed that the solution state is \a z n = 1} with no loss of generality. The factor 
of N 2 ensures that the energy scale of all terms in the Hamiltonian are of order ~ N 2 , which 
can be executed using the quantum bus methods as described above. The bosonic qubits are 
prepared in the state \X) = Yln=i It^> 7l)) n an< ^ ev °l ve d in time by applying H. The system 
then executes Rabi oscillations between the initial state \X) and the solution state \ANS). 
The time required for reaching the solution state can be estimated from the the period of 
the Rabi oscillations. For a system undergoing Rabi oscillations of a form (S z )/N = sin 2 ut, 
the frequency can be estimated by to = \J d ^ S J^' Z ^-/2N. Evaluating the second derivative 
in the Heisenberg picture gives d ^jgr^ = — (Wc [Hg, S*/N]\) = 2N 2 /2 M , corresponding to 
a evolution time of t ~ \Jl M /N '. The bosonic version of the algorithm has the same square 
root scaling with the number of sites, but with a further speedup of N, resulting from the 
ultrafast gates made possible by the use of bosonic qubits. A numerical calculation for a 
simple two site case is shown in Figure Eb, which clearly shows the factor of N improvement 
in speed of the Grover algorithm. 

In conclusion, we have found that using qubits made of bosons occupying two level 
systems connected by a quantum bus can realize two qubit gates at a boosted energy scale 
of factor iV 2 with no penalty on decoherence times. This allows for gate times to be reduced 
by a factor N, resulting in many more quantum operations in a given decoherence time. 
The most likely realization of such two level bosonic systems are in BECs which can be 
manipulated using optical transitions. Our scheme differs from previous proposals of using 
3ECs as qubits in the way the states are encoded in the bosons. Unlike previous proposals 
]jj where fragile Schrodinger cat-like states are used, the encoding presented here is robust 
against decoherence. For unitary rotations, there is a simple mapping connecting standard 
qubits to the bosonic qubit states, giving an equivalent operation. We leave applications 
including non-unitary operations such as quantum teleportation as future work. 
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